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The effect of the electron-electron cusp on the convergence of configuration interaction (CI) wave 
functions is examined. By analogy with the pseudopotential approach for electron-ion interactions, 
an effective electron-electron interaction is developed which closely reproduces the scattering of 
the Coulomb interaction but is smooth and finite at zero electron-electron separation. The exact 
many-electron wave function for this smooth effective interaction has no cusp at zero electron- 
electron separation. We perform CI and quantum Monte Carlo calculations for He and Be atoms, 
both with the Coulomb electron-electron interaction and with the smooth effective electron-electron 
interaction. We find that convergence of the CI expansion of the wave function for the smooth 
electron-electron interaction is not significantly improved compared with that for the divergent 
Coulomb interaction for energy differences on the order of 1 mHartree. This shows that, contrary to 
popular belief, description of the electron-electron cusp is not a limiting factor, to within chemical 
accuracy, for CI calculations. 



I. INTRODUCTION 

The primary difficulty in solving the Schrodinger equa- 
tion for many-electron systems arises from the presence 
of the electron-electron interaction, which results in a 
problem non-separable in the electron coordinates. Con- 
figuration interaction (CI) calculations use an expansion 
of the many-body wave function in configuration state 
functions (CSF) given by spin coupled sums of determi- 
nants built from single-particle orbitals. One of the cen- 
tral difficulties for the convergence of the CI expansion 
is thought to arise from the behaviour of the Coulomb 
interaction at small electron-electron separations [^^^Oj. 
The many-body wave function must have an appropriate 
cusp at electron coalescences so that the infinite Coulomb 
\ interaction term is exactly cancelled by an opposite di- 
vergence of the kinetic energy . 

It is often stated in the literature that failure of the 
CI expansion to reproduce the correct electron cusp, i.e. 
the short-range part of the Coulomb hole, leads to a slow 
convergence in the energy with respect to the number 
' of CSF's §. While this has been shown [||] in the 
asymptotic regime (i.e. very large angular momentum I 
and energy errors of the order of 10^^ Hartree), it is not 
clear that it applies in the practical regime (mHartree 
accuracy) for multi-electron systems. 

Explicit inclusion of terms in the interelectronic dis- 
tance (ri2 methods) [p^ji) or the use of correlated (gem- 
inal) Gaussian basis functions significantly improve 
the convergence of a CI calculation. For the helium atom, 
the inclusion of ri2 terms Q yields an expansion in partial 
waves converging asymtotically as or better, instead 
of a slow , where I is the maximum angular momentum 
in the expansion. Baker et al. also showed that using 
basis functions with the same analytic structure as the 
true wave function (depending on powers and logarithms 
of ri2) leads to a reduced expansion length. The suc- 
cess of methods explicitly including ri2 terms has been 



attributed to their correct description of the short-range 
part of the interaction. 

A completely different point of view, taken by Gilbert 
in 1963 is that the difficulty of CI wave functions in 
describing short-range effects of the cusp are energetically 
unimportant and that the benefit of the explicitly corre- 
lated approach is that the entire Coulomb hole, with a 
size on the order of the atomic orbitals, is much easier to 
describe in terms of interelectronic coordinates. Gilbert's 
suggestion appears to have received little attention and 
the idea that smoothing out the Coulomb potential and 
removing the cusp altogether will automatically lead to 
an improved CI convergence has been stated in the liter- 
ature |l6|j9t|. However, no systematic studies of the effect 
of the electron cusp on CI energy convergence at the level 
of interest for chemical accuracy (1 mHartree for total 
energies) have been presented. 

In this paper, we explore whether a significant im- 
provement in the convergence of the CI expansion is 
obtained when the true Hamiltonian is replaced by one 
which has a smooth and finite electron-electron interac- 
tion (pseudo-interaction) at short distances and, conse- 
quently, smooth eigenfunctions with no cusps. 

Surprisingly, in calculations for the He and Be atoms, 
we find that removing the divergence of the Coulomb 
potential at electron-electron coalescences does not sig- 
nificantly reduce the length of the CI expansion, thereby 
substantiating the argument put forward by Gilbert p^ ] 
that the slow convergence should be attributed to dif- 
ficulties in describing, not the interelectron cusp, but 
the intermediate-range electronic correlation. The poor 
treatment of the interelectron cusp in the CI expansion is 
not a limiting factor at the mHartree level of accuracy on 
the total energy. To further confirm this result, we calcu- 
late the energy within quantum Monte Carlo (QMC) for 
the true electron-electron interaction with accurate wave 
functions which explicitly contain interelectronic coordi- 
nates, and then deliberately remove the short-range cor- 
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relation term related to the cusp, keeping other parts 
of the wave function unchanged. We find that the re- 
moval of the cusp has a very small eS^ect on the vari- 
ational energy. Finally, we demonstrate that the con- 
vergence behavior of the CI expansion as a function of 
angular momentum in the range practical for many ap- 
plications (Z < 6) is largely unaffected by the smoothing 
of the electron cusps. Therefore, the angular momentum 
convergence numerically observed by other authors is 
related to problems associated with intermediate range 
correlations, not the cusp region as such. 

In Sec. 0, we describe the pseudopotential used for 
the electron-electron interaction and briefly present the 
QMC and CI methods. In Sec. |ll| we report the QMC 
and CI energies for He and Be atoms using the true and 
pseudo electron interaction. In Appendices ^ and 
we discuss the numerical generation and accuracy of the 
electron pseudo-interaction, and in Appendix |^ we define 
a short-range correlation term which can be associated 
with it. In Appendix^ we briefly discuss the evaluation 
of two-electron integrals for the pseudo-interaction. 



II. COMPUTATIONAL METHODS 
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FIG. 1. Norm-conserving electron pseudo-interaction 
(solid line) and smooth interaction Krf (dashed line), each 
with a cut-off radius rc = 0.25 a.u. (indicated by a vertical 
dotted line), and the true 1/r interaction (dash-dot line), as 
functions of the electron separation r. 



A. Electron-electron Pseudo-interaction 

When modifying the short-distance electron-electron 
interaction, we wish to keep the intermediate- and long- 
range correlation properties of the many-body eigen- 
states unchanged. This will allow us to isolate systemati- 
cally the effects of the short-range cusp-like behaviour in 
the wave function. To achieve this, our new (smooth) in- 
teraction should ideally have electron-electron scattering 
properties identical to the true interaction for all ener- 
gies. If this were the case, all many-electron eigenvalues 
would be identical for both interactions. While this ideal 
cannot be exactly realized, we expect that a potential 
with the same scattering properties as the true interac- 
tion over a sufficiently wide range of scattering energies 
will give many-electron eigenvalues very close to the true 
ones. 

This goal for electron-electron scattering is analogous 
to the problem of electron-ion scattering addressed in the 
well-established pseudopotential method . In that ap- 
proach, the aim is to produce a smooth electron-ion po- 
tential which reproduces the single-particle eigenvalues 
of the true system over a range of energies, and has the 
same single-particle wave functions outside some suitably 
chosen cut-off radius r,,. Such pseudopotentials are gen- 
erated for the isolated atom but may then be used with a 
high level of accuracy in molecules and solids. In solids, 
the smoothness of the potential greatly improves the con- 
vergence of plane-wave expansions for single-particle or- 
bitals |l8|. 

Following the success of norm-conserving pseudopo- 
tential methods in the context of the electron-ion scat- 



tering problem, we generate a norm-conserving pseudo- 
interaction for the two-electron scattering problem, 
which will replace the true e^/r electron-electron interac- 
tion inside a cut-off radius Tc- Since the electron-electron 
interaction is repulsive and has no bound states for an 
isolated pair of electrons, we use, with some minor mod- 
ifications, the generalized norm-conserving method de- 
vised by Hamann to generate pseudopotentials for 
electron-ion interactions at unbound, scattering state en- 
ergies. The numerical details of this approach are pre- 
sented in Appendix ^ 

In Fig. 0, we show a pseudo-interaction potential V-ps 
generated by this approach, along with the true e^/r in- 
teraction and a smooth but not norm-conserving inter- 
action, VeTi(r) = e^eYi{r/rc)/r, as used by Savin and 
coworkers in a different context. In Appendix |b[ the 
accuracy of such scattering potentials in reproducing the 
exact eigenstate energies is analyzed qualitatively using 
first-order perturbation theory. In Appendix |c|, we de- 
fine a short-range inter-electronic correlation term, which 
can be directly related to the replacement of the divergent 
Coulomb interaction with a smooth norm-conserving in- 
teraction. 



B. Quantum Monte Carlo Calculations 

We employ QMC methods to determine accurate en- 
ergies both for the Coulomb and the pseudo-interaction. 
The wave functions used in QMC are of the form given 
in Ref. [2l| and explicitly depend on the inter-electronic 
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coordinates. They include parameterized terms to de- 
scribe electron-nucleus, electron-electron and electron- 
electron-nucleus correlation, as follows: 

(Y.dnBlBl}j3{r,,r,,n,) . (1) 

DT and are the Slater determinants of single particle 
orbitals for the spin-up and down electrons, respectively, 
and J is a Jastrow factor of the form, 

J(r^,rj,ry) = J| exp(^i) J| exp(By) J| exp(Py), (2) 

i ij ij 

where A is an electron-nuclear term (which could be 
omitted if a sufficiently large single particle basis were 
used in constructing the determinantal part), B is an 
electron-electron term incorporating the cusp at small 
Vij, and P is a smooth function of the electron coordi- 
nates, Ti and discussed in detail in Ref. For 
the pseudo-interaction calculations, the B term is omit- 
ted. The parameters in the determinant and Jastrow 
factor are optimized within QMC using the variance min- 
imisation method and the accuracy is tested at the 
variational level. The wave function is then used in dif- 
fusion Monte Carlo (DMC), which produces the best en- 
ergy within the fixed-node approximation (i.e. the lowest- 
energy state with the same nodes as the trial wave func- 
tion). 

For a given interaction, this method yields the ex- 
act eigenvalue (within statistical sampling noise) for the 
ground state of the He isoelectronic series since the sin- 
glet ground state wave function of a two-electron system 
is nodeless. For He, only one determinant is used in the 
wave function. For the Be atom, the 2s and 2p configura- 
tions are included in the determinantal part of the wave 
function and DMC gives the best energy subject to the 
fixed-node condition. To isolate the effect of the pseudo- 
interaction from the fixed-node approximation, we keep 
the nodal structure of the wave function fixed while mod- 
ifying the electron-electron interaction. 

C. Configuration Interaction Calculations 



For use in the CI approach, the pseudo-interaction Vps 
is represented (in Hartree atomic units) as a sum of the 
smooth function Verf plus no Gaussian terms, as follows: 

Fp3(r) = ^I!MIA + f^c,eM-d^r') ■ (3) 

i=l 

This reproduces the l/r tail of the Coulomb interaction 
and the smooth, finite form of the pseudo-interaction at 
short distances. The parameters Ci and di are found by 
least-squares fitting (for radii r < 3rc) to the numerical 
potential, generated as described in Appendix 0. For 
the CI results presented below, the potential was fitted 
using riQ = 9. QMC calculations of the He atom energy 
with the numerical potential and with the fitted potential 
show a negligible difference (less than 0.02 mHartree for 
all values of Tc < 1 a.u.). 

Evaluation of the pseudo-interaction Vps (r) can be im- 
plemented in codes for Gaussian integral calculation since 
it is possible to compute the two-body matrix elements 
both of exp(— d r^) and eii{r/rc)/r in the CI basis. Ma- 
trix elements of the interaction evi{r/rc)/r for all s-type 
basis functions have already been given in Ref. For 
both types of interaction, a simple modification to the 
integral codes is possible if the standard two-electron in- 
tegrals are calculated using the method of Rys quadra- 
ture Details of two-body integral evaluation with 
the modified form of the electron interaction are given in 
Appendix 

III. He AND Be ATOMIC ENERGIES WITH 
PSEUDO-INTERACTIONS 

To test the impact of the interelectron cusp on the en- 
ergy, we calculate the total energies of He and Be atoms, 
replacing the true 1 /r electron interaction with the norm- 
conserving pseudo- interaction Vps(r), generated for var- 
ious values of the cut-off radius re- For comparison, 
we also present results for the smooth but not norm- 
conserving potential VaTiir) = evl(r / r / r . The results 
in this section are obtained from QMC and CI calcula- 
tions. 



The CI calculations are performed using the Monte 
Carlo configuration interaction (MCCI) method p3| , p4t . 
In a MCCI calculation, a coefficient threshold Cmin is de- 
fined, configurations are randomly generated, and, after 
a variational calculation, coefficients of magnitude be- 
low Cmin are rejected. This procedure is repeated until 
a desired degree of convergence in the CI calculation is 
achieved. A full CI (FCI) calculation can be obtained for 
Cmin = 0. The fact that the MCCI procedure does not 
rely on a pre-selection of CSF's makes it well-suited for 
our study: the dynamical selection of configurations al- 
lows for rejection of CSF's which may no longer be needed 
with the smooth electron interaction. 



A. Ground State of tlie He atom 

In Fig. |[ we show the singlet ground state energies of 
the He atom calculated with a norm-conserving electron- 
electron interaction Vps generated as in Sec. for various 
values of the cut-off radius Tc- The values obtained from 
QMC (exact) calculations are shown. We also plot the 
energies calculated using the V^i [r) interaction. The hor- 
izontal line indicates the exact energy Eq — —2.903724 
Hartree for the true Coulomb interaction |^ . 

The calculated QMC energies for Vps and Krf in 
Fig. 1^ are fitted with polynomial expansions of the form: 
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FIG. 2. He atom ground state energy i5 as a function 
of cut-off radius rc, obtained from DMC calculations using 
various forms of the electron-electron interaction: (a) the 
norm-conserving pseudo-interaction Vps (solid line); (b) the 
smooth interaction 14rf (dashed line). The curves are poly- 
nomial fits to the calculated DMC values, as discussed in the 
text. The horizontal line shows the exact value for the He 
ground state. Also shown are variational energies, with the 
true Coulomb interaction, upon removal of the electron cusp 
(see Appendix ^) from the wave function used in the DMC 
calculations (dotted line). 



E{rc) = Eq + arl + &rj? -t- cr^ -f dr^ -f er^. As discussed 
in Appendix any potential which differs from the ex- 
act 1/r interaction would be expected, from first order 
perturbation theory, to converge at least quadratically in 
Tc to the exact energy Eq for small Tc- Since Krf is not 
explicitly constructed to eliminate first order errors in 
yc-[i{f) — l/?", quadratic convergence is seen for Vc-d and 
the fitting parameters are: a = —0.3279, b = —0.3839, 
c = 0.6758, d = -0.2622, and e = 0. On the other 
hand, the no rm-c onserving pseudo-interaction construc- 
tion of Sec. [I A explicitly ensures that the first-order 
perturbation correction to the energy is zero. The He 
ground state energy as a function of using Vps can 
indeed be fitted without the quadratic term, and even 
the cubic term proves to be small, the polynomial coef- 
ficients being: a = 0, b = 0.0188, c = -0.1023, d = 0, 
and e = 0.0280. For He, the energies obtained using 
Vps(r) are within 1 mHartree of the energy computed 
with the true Coulomb interaction for values of Tc < 0.4 
a.u. Since the typical electron-electron distance in the 
He ground state is of the order of 1 a.u., these results in- 
dicate that the first peak of the electron pair-distribution 
function provides a good measure of how large a cut-off 
radius Tc may be used to obtain accurate total energies 
with Vps{r): rc should be approximately half the radius 
at which this first peak occurs. The same criterion for 
the largest allowed Vc value was found to hold for other 
systems (not shown), viz. the He triplet state, the He 
isoelectronic series, and multi-electron atoms. In a qual- 



itative way, this condition is similar to the criterion used 
in the generation of electron-ion pseudopotentials, where 
the cut-off radius is chosen to be inside the peak of the 
valence electron density. 

Also shown in Fig. || are the variational energies (with 
the true Coulomb electron-electron interaction) obtained 
for the He ground state by removing the electron cusps 
from the accurate trial wave function which was deter- 
mined by variance minimiza tion for use in the DMC cal- 
culations, as specified in Eq. |Cl| . This trial wave function 
includes considerable variational freedom in the interelec- 
tronic terms and obtains 100% of the correlation energy. 
As can be seen in Fig. ^, removal of the cusp out to radii 
of 0.5 a.u. causes little degradation of the quality of the 
energy obtained with this wave function. Thus the suc- 
cess of interelectronic terms in obtaining good variational 
energies is not primarily due to their description of the 
short-range cusp behaviour. 

This trend is further borne out when we look at 
the convergence of partial CI expansions within a given 
single-particle basis set for the MCCI method. The CI 
calculations were performed with the correlation consis- 
tent polarized valence double-zeta (cc-pVDZ), triple-zeta 
(cc-pVTZ), quadruple-zeta (cc-pVQZ) and quintuple- 
zeta (cc-pV5Z) basis sets of Dunning with the elec- 
tron interaction Vps for various values of rc- All four 
basis sets show similar behaviour and we present only 
the cc-pVQZ results in detail. Table | shows the correla- 
tion energy obtained as a function of cut-off radius for 
various values of the expansion coefficient threshold Cmin- 
The full CI result is obtained when Cmin — 0.0. For larger 
values of the cut-off radius (upper rows of table) , the 
total correlation energy is smaller, as the variation of the 
interaction Vps becomes smaller around the typical sepa- 
ration distance of the electrons. However, the fraction of 
the correlation energy which is recovered at a given CI 
threshold level Cmin depends little on the value of re- 
Fig. |3| shows the number of CSF's in the CI expan- 
sion as a function of cut-off radius rc for various values 
of the expansion coefficient threshold c^in- For values of 
Cmin > 3 X 10~^, One is quickly reduced to a few CSF's 
and finally to the Hartree-Fock solution. As with the cor- 
relation energies, for a given threshold Cmin, the number 
of CSF's shows relatively little change for values of rc be- 
tween and 0.95. Again, it is worth bearing in mind that 
the typical electron-electron separation is 1 a.u. and that 
a cut-off radius rc = 0.95 represents a very smooth inter- 
action. Thus, removal of electron cusps does not result 
in a significant reduction in the number of contributing 
configurations in the CI expansion. 

We next demonstrate that removal of the electron cusp 
does not significantly change the convergence properties 
of CI expansions as the highest angular momentum quan- 
tum number is increased in the single particle basis sets. 
In Fig. ^, FCI energies for He are displayed against in- 
creasing basis set size in the correlation consistent series 
from Dunning, and hence increasing angular momentum 
number in the basis. The convergence in the FCI energy 
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TABLE I. Calculated correlation energy for the He atom with a cc-pVQZ basis set, as a function of cut-off radius 
and the CI coefficient threshold Cmin- The energy units are mHartree. The right-hand column (HF— QMC) shows the full 
correlation energy (i.e. the difference between the Hartree-Fock energy, calculated with this basis, and QMC (exact) energy). 
The percentages indicated are the fractions of the full CI correlation energy for each value of rc. 



(a.u.) 


3 X 10"^ 


2 X 10"^ 


1 X 10"^ 


1 X IQ-^ 


1 X 10"" 


0.0 


HF - QMC 


0.95 


2.15 


15.10 


12.98 


18.72 


18.96 


18.96 


19.56 




11.01% 


77.21% 


66.36% 


95.69% 


96.94% 


96.95% 


100.00% 


0.80 


6.09 


18.45 


16.93 


25.11 


25.16 


25.16 


25.76 




23.66% 


71.64% 


65.74% 


97.48% 


97.67% 


97.67% 


100.00% 


0.50 


8.64 


27.23 


32.79 


36.32 


36.32 


36.32 


37.18 




23.24% 


73.24% 


88.19% 


97.69% 


97.69% 


97.69% 


100.00% 


0.20 


14.07 


28.93 


35.27 


40.75 


40.75 


40.75 


41.93 




33.57% 


69.00% 


84.13% 


97.20% 


97.20% 


97.20% 


100.00% 


0.10 


9.37 


28.98 


35.35 


40.89 


40.89 


40.89 


42.18 




22.22% 


68.71% 


83.80% 


96.94% 


96.94% 


96.94% 


100.00% 


0.00 


27.44 


30.61 


36.09 


40.90 


40.90 


40.90 


42.21 




65.01% 


72.52% 


85.50% 


96.90% 


96.90% 


96.90% 


100.00% 



60 |— r 




(a.u.) 

FIG. 3. The number of configuration state functions 
(CSF) in the CI expansion for the He atom with the cc-pVQZ 
basis set as a function of Tc- Results are presented for coeffi- 
cient thresholds Cmin of 0.03, 0.02, 0.01, 0.001 and 0.0001. 
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FIG. 4. Full CI energies E for the ground state He atom 
as a function of L, the maximum angular momentum in the 
basis set used, calculated using the electron interaction Vps 
with values of the cut-off radius rc from to 0.8. 



with increasing basis set size is compared for electron in- 
teraction Vps, with values of the cut-ofF radius Tc from 
to 0.8. Fitting these energies to a function similar to the 
asymptotic form derived in Ref. H, 



E 



a b 



(4) 



where L is the maximum angular momentum in the basis, 
gives best values of a ~ 0.038 and h « —0.022 for each 
value of Tc < 0.2 and gives a ~ 0.02 and b k, —0.01 for 
Tc > 0.5. Thus, the predominant convergence of the 
energy is present for all values of rc investigated here. 
Although the asymptotic form in Eq. ^ fits the results 



well, the parameters a and b do not have the asymptotic 
values, a = 0.025 and b = 0.008, found for the angular 
momentum convergence associated with the short-range 
cusp 1^. Two factors cause this discrepancy: (i) not only 
higher angular momentum functions are added as L is in- 
creased in the Dunning basis sets; (ii) more importantly, 
we are not near the asymptotic regime, where expansion 
of the short-range cusp would be the dominant energy 
correction. 
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FIG. 5. Be atom ground state energy i5 as a function of 
cut-off radius Vc, obtained from (a) DMC calculations using 
the norm-conserving pseudo-interaction Vps (solid line); (b) 
variational calculations using the true Coulomb interaction 
upon removal of the electron cusp from the trial wave function 
(see Appendix ^ used in the DMC calculation (dashed line) . 



TABLE II. Calculated correlation energy, as a function of 
cut-off radius rc and CI coefficient threshold Cmin, for the Be 
atom with the TZP basis set. All energies are in mHartree. 
The right-hand column shows the QMC correlation energy 
(i.e. the difference between the Hartree-Fock energy for this 
basis and the QMC energy) for each value of the cut-off radius 
Vc- The percentages indicated are the fractions of the QMC 
correlation energy for each value of r^. 
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TABLE III. Calculated correlation energy for the Be atom 
with the ANO basis set. The notation is the same as in Ta- 
ble 0. 



B. Ground State of the Be atom 

The qualitative nature of the CI energy convergence 
seen in the previous section is not special to the He 
atom. Fig. [s] shows the DMC calculated ground state 
energies of Be for various cut-off radii rc- Because of the 
much smaller radius of the Is orbital in Be, compared to 
He, the energetically appropriate cut-off radius Vc for the 
pseudo-interaction is substantially smaller. Also shown 
in Fig. H are the variational energies obtained for the Be 
ground state by removing the electron cusps (as defined 
in Eq. CI) from the accurate trial wave function which 
was determined by variance minimization for use in the 
DMC calculations. This trial wave function includes in- 
terelectronic terms and obtains 99.5% of the correlation 
energy. Removal of the cusp causes significant degrada- 
tion of the quality of the wave function only for rc > 0.3 
atomic units, which is comparable to the typical sepa- 
ration of the Is electrons. Thus, in Be, as in He, the 
success of interelectronic terms in obtaining good varia- 
tional energies is not primarily due to their description 
of the short-range cusp behaviour. 

CI calculations using basis sets with s- and p- functions 
(the TZP basis of Ahlrichs et al. j28|), s-, p- and d- 
functions (6-311g* ||2^]), s-, p-, d- and f-functions (the 
ANO basis of Widmark et al. 30 1) and s-, p-, d-, f- and 
g-functions (i.e. the ANO basis, augmented with a g- 
function with exponent 1.2) were carried out. The values 
of the CI coefficient threshold Cmin were taken to be 10"^, 
IQ-'^ and 10"'' for all basis sets. Calculations were car- 
ried out at cut-off radii Tc of (i.e. Coulomb interaction), 
0.1 and 0.2 a.u. 
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In Tables ttl and [II we present the correlation energy 



for Be, obtained with the TZP and ANO basis sets, re- 



spectively. Fig. ^ shows the number of CSF's as a func- 
tion of the cut-off radius rc for various values of the CI 
coefficient threshold Cmin in the MCCI approach. We see 
that, especially at the higher levels of accuracy (i.e. small 
Cmin), neither the fraction of the correlation energy ob- 
tained nor the number of CSF's used depend strongly on 
the cut-off radius rc- As we saw in the He atom, remov- 
ing the divergence of the Coulomb electron interaction 
does not substantially change the CI energy convergence 
at this level of accuracy. 



IV. CONCLUSIONS 

In this paper we have directly explored the effect 
of the electron cusp on the convergence of the en- 
ergy for CI wave functions. To do so, we introduced 
a fictitious electron-electron interaction which, unlike 
the true Coulomb potential, does not diverge at elec- 
tron coalescences and therefore has many-body eigen- 
functions which are smooth there. The smooth poten- 
tial is obtained by a method analogous to the well- 
established treatment of electron-ion interactions with 
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range cusp. Rather, the improvement in convergence at 
the mHartree level must be understood in terms of a 
better description of medium-range correlations when ex- 
pressed in interelectronic coordinates. This fact must be 
considered when developing theoretical treatments of the 
electronic correlation problem. 
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FIG. 6. The number of configuration state functions 
(CSF) in the CI expansion for the Be atom with TZP and 
ANO basis sets. Results are presented for coefficient thresh- 
olds Cmin of 0.01, 0.001 and 0.0001 at cut-off radu Tc of 0.0, 
0.1 and 0.2 a.u. 



smooth, norm-conserving pseudopotentials. These give 
scattering properties very similar to the true potential 
and improve the convergence of plane-wave expansions 
for single-particle wave functions. By analogy, one might 
expect that a similar treatment of the electron-electron 
interaction would improve the convergence of CI expan- 
sions. 

For the He and Be atoms, we have compared the CI 
expansions of wave functions for the true, Coulomb inter- 
action and for smooth pseudo-interactions between the 
electrons. We have also compared the variational en- 
ergies of highly accurate wave functions containing vyi 
terms, where the cusp condition is exactly satisfied, with 
the energy for the same wave functions, from which the 
cusp has been deliberately removed out to a given radius. 

The main result of our study is that a description of 
the electron cusp as such is not a limiting factor in cal- 
culating correlation effects with configuration interaction 
methods at the mHartree level of accuracy. We are led 
to this conclusion since the replacement of the diver- 
gent Coulomb interaction with a finite interaction leaves 
the convergence properties largely unchanged. The slow 
convergence of CI expansions in this energy range must 
be attributed to medium-range correlations, which are 
present for both types of electron interaction. The results 
obtained by explicitly removing the cusp from accurate 
variational wave functions support this conclusion, show- 
ing that the cusp of the wave function is not energetically 
important. This is found, even when cusp-like behaviour 
is removed out to a distance of half the typical interelec- 
tronic separation. 

When describing correlation effects with explicit r\2 
methods, it is a mistaken notion to assign the improve- 
ment in convergence to their ability to describe the short- 



APPENDIX A: GENERATION OF ELECTRON 
PSEUDO-INTERACTIONS 

The procedure for constructing the electron pseudo- 
interaction Vps follows Hamann's construction of gen- 
eralised norm-conserving pseudopotentials for electron- 
ion scattering, with two modifications: (i) Only two 
particles are involved in the generation of the electron- 
electron pseudo-interaction, so that we do not need to in- 
clude screening potentials; (ii) because we are solving the 
scattering problem in the centre of mass frame of the elec- 
trons and use the relative coordinate r of the electrons, 
the mass of the electron is replaced by the reduced mass 
for the two-electron system. In Hartree atomic units, the 
interaction potential is V(r) = \/r and the reduced mass 
is 1/2. The radial Schrodinger equation for the relative 
motion of two electrons is then 



(Pui ^ + 



+ {V ^ei)\ui = Q, 



(Al) 



where ui{r) is r times the wave function in the relative 
coordinates, with angular momentum / and energy ti (the 
kinetic energy of the separated electrons in the centre of 
mass frame). 

In the generation of a pseudopotential, the true po- 
tential and pseudopotential become identical outside a 
suitably chosen cut-off radius Vc and are constructed to 
have identical scattering phase shifts and energy deriva- 
tive of the phase shift at a reference energy ti . A smooth 
pseudo wave function u-psir) is generated numerically to 
ensure "norm-conservation", i.e. Ups{Roo) — ui{Roo) and 



Wps(r)^dr 



(A2) 



for all Roc >> Tc- In practice, we choose Rao ^ 2.5rc 



as in Ref. 1^1 f 



Given 



a smooth, finite pseudopo- 



tential Vps is generated by inversion of the Schrodinger 



Eq. Al for r < Roo, with ui and V substituted by Ups and 



Vps, respectively. We follow the steps exactly as given in 
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Sec. II of Ref. |T^], except where the reduced mass en- 
ters in the final inversion, sUghtly modifying Eq. 12 of 
Hamann's work. 

In the norm-conserving pseudopotcntial approach, dif- 
ferent scattering potentials are often generated for differ- 
ent values of the angular momentum I. The full many- 
body wave function is anti-symmetric with respect to ex- 
change of any two electrons. Thus, if the two electrons 
have parallel spin, the orbital part of the wave function 
must have odd parity and the allowed angular momentum 
I must be odd. Thus, the dominant parallel-spin scat- 
tering has 1 = 1. For singlet (anti-parallel) scattering, 
the orbital part is symmetric and / must be even, with 
dominant allowed angular momentum 1 = 0. Following 
the usual pseudopotcntial approach, we might then ex- 
pect to generate a parallel-spin scattering potential using 
I = 1 and an anti-parallel-spin scattering potential using 
1 = 0. In practice, we find that the I = potential gives 
good scattering for both angular momenta and the use 
of separate parallel- and anti-parallel-spin potentials is 
unnecessary. 

Since e; is positive, the wave function will oscillate for 
large r, and if a node occurs for r < Roo, the inver- 
sion of the Schrodinger equation will fail. Thus, the joint 
requirements, that Tc be sufficiently large to produce a 
smooth pseudo-interaction and that Rao ^ 2.5rc, effec- 
tively limit the maximum reference energy e; at which 
the pseudo-interaction construction is possible. 

We find that the pseudo- interaction is transferable (i.e. 
has very similar scattering strength to the true interac- 
tion) for energies over a wide range about the reference 
energy e; , as found previously for electron-ion pseudopo- 
tentials. In particular, we find that transferability is best 
for energies less than the reference energy at which the 
pseudopotcntial was generated. This would suggest that 
the best reference energy to use is the largest energy 
compatible with the procedure for inversion of the radial 
Schrodinger equation, a choice confirmed by calculations 
of many-body eigenstates using pseudo-interaction po- 
tentials generated with different values of the reference 
energy. 

APPENDIX B: FIRST ORDER PERTURBATION 
ANALYSIS OF PSEUDO-INTERACTIONS 

In this Appendix, we consider the first-order error 
in the energy expected for a typical, arbitrary smooth 
electron-electron potential which is equal to the exact in- 
teraction 1/r outside some (adjustable) cut-off radius rc- 
We will assume that the potential is (at least approxi- 
mately) of the form 

VrAr) = -U{r/r,) (Bl) 
rc 

where U{r/rc) is some universal function independent 
of rc. Clearly both the true interaction 1/r and the 



smooth interaction T4rf are of this form exactly. Vps 
is also approximately of this form. Thus, the difference 
Ay between the interaction Vr^ (r) and the true electron- 
electron interaction 1 /r is of the same form: 

AV{r) = -AU{r/rc), (B2) 

rc 

where AU{p) is zero for p> 1. 

Replacing the true electron-electron interaction with 
Vr^ (r) then gives rise to a first-order energy error: 

6e^^'> = (*|At/|*) 

where \'^) is the eigenstate of the many-electron system. 
If P(r) is the electron pair distribution function (aver- 
aged at the pair separation distance r), the energy error 
may be rewritten as: 

f°° 1 
5e^^^ = / P{r) — AU{r/rc)A7Tr^dr 
Jo rc 

= rl ( P{rcp)AU{p)4np''dp. (B4) 
Jq 

If rc is sufficiently small and the pair-distribution func- 
tion P{r) does not tend to zero at r = 0, then P{rcp) may 
be replaced by P(0) in the integral and 5t^^^ ~ KP(0)r1 
for small rc, where K \s a, constant for the type of inter- 
action, defined by 

K = [ AU{p)4Trp^dp. (B5) 
Jq 

The construction of Vps ensures that the energy error 
is identically zero for the two-electron scattering at the 
reference energy, and therefore K = for Vps. However, 
for an arbitrary smooth interaction (such as Vert), K will 
not be zero. 



APPENDIX C: SHORT-RANGE CORRELATION 
TERM OBTAINED FROM NORM-CONSERVING 
POTENTIAL 

In the generation of the norm-conserving potential dis- 
cussed in Appendix we obtained a pseudo wave func- 
tion Ups(r), which is a solution to the radial Schrodinger 
equation at the scattering energy e; for the pseudo- 
interaction, and the wave function u;(r), which is the 
corresponding solution for the true Coulomb interaction. 
The ratio J^r = ui/ups of these two functions has the 
correct electron cusp (i.e. logarithmic derivative of —1/4 
or —1/2 for parallel or anti-parallel spins, respectively, at 
r = 0) and tends rapidly to 1 for r > r^. This function J^r 
may be said to contain the short-range cusp behaviour 
of the Coulomb scattering solution, since it represents 
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the difference between tlie scattering wave function for 
the divergent potential and that for the finite pseudo- 
interaction. 

Moreover, if ^'(ri, . . . , tat) is a many-electron wave 
function (e.g. one which has inter-electronic coordinates 
rij explicitly included) with the correct electron cusps, 
the function 



*s = *(ri,...,rAr)/J|j^^(ry) 



(CI) 



i<j 



is smooth at electron coalescences but has very similar 
electron correlations for > re- The wave function 
has the electron cusp behaviour of the original wave 
function removed inside the cut-off radius Tc, but is 
otherwise the same as ^E*. Thus, the energy difference 
(^'s|-ff|^'s) — ('J'l^l^) is a good measure of the energy 
cost of removing the electron cusp behaviour for r < 
from a correlated wave function, while keeping other 
aspects of the wave function (e.g. any benefits there may 
be to the inclusion of inter-electronic coordinates in the 
form of the trial function) unchanged. 



The Rys scheme represents the l/ri2 term in Gaussian 
integral form, to obtain: 

1 2 f°° 

((?!)i0j| — \(l)k(f>i) = —1= / du{(f>i(l)j\exp{-u'^rl2)\4ik4'i)- 
ri2 V"' Jo 

(Dl) 

This integral is then transformed to an integral over the 
variable t, where = u^/ip + u^) and 1/ p = !/(«; -|- 
aj) + l/{ak + ai): 



h4>j\ — \4>k(t>i 

T\2 



[ PL(t)exp{~Xt^) dt , (D2) 
Jo 



where Pl is a polynomial and X = p\rA — fsp. Here, 
rA and rs are the usual weighted- average centers of the 
Gaussian orbital pairs and {k,l), respectively, as 

given in Ref. . The integration over t is then evaluated 
exactly by an n-point quadrature formula: 



) = ^PL(ta)W„ 



APPENDIX D: EVALUATION OF 
TWO-ELECTRON GAUSSIAN INTEGRALS 
WITH PSEUDO-INTERACTIONS 

In this Appendix, we give details of the evaluation of 
two-electron Gaussian basis set integrals for the electron 
interactions, V^rf and Vps, and their incorporation into 
existing electron integral programs. 

For the interaction Vgrf, the necessary modification of 
two-electron integrals for all s-type basis functions have 
been given by Savin |2^]. We do not use Savin's ap- 
proach, but rather directly modify the Rys quadrature 
approach of Dupuis, Rys, and King |2^], to obtain ma- 
trix elements for the interaction Veii in place of those for 
l/?'i2- To evaluate two-electron integrals for the pseudo- 
interaction V^s, we use the form given in Eq. |3[ This in- 
volves the evaluation of integrals with the interaction Vcrf 
and a sum of Gaussian terms. In codes (e.g. the ARGOS 
code iQ) which use the Rys quadrature approach for 
the evaluation of two-electron integrals with the Coulomb 
interaction, these Gaussian terms are already evaluated 
in an intermediate stage and may be used in the present 
context. 

Let us define two-electron integrals of the Coulomb 
interaction as 



ri2 



4>i{ri)4>j{ri) (I>k{r2)(f>iir2)dridr2 

ri2 



for Gaussian basis functions, 



,(r) = exp(-a,|r - R,\^) Y[{x- X,)"^ 



where ta{X) and Wa{X) are determined as in Ref. p5|| . 

We observe that the interaction Vj^ rf m ay be written 
in exactly the same form as in Eq. Dl, replacing the 
upper limit of integration (u = oo) with u ~ l/rc. This 
changes the upper limit of integration in the variable t 
from 1 to = 1/ yl 



pr'^ in Eq. D2. To use directly 



the Rys quadrature approach, we further transform to 
the variable t' — t/tc, to obtain 



((Pi(Pj\ \<Pk<Pl) = I. 

ri2 

which can be evaluated as. 



; [ PLitct')exp{~Xtlt'^) dt' 
Jo 



erf(ri2/rc) 

Djl \(Pk<Pi 

ri2 



where the quadrature points t'^ and weights are 
found by defining X' = Xt^, t'^ = tcta(X'), and = 

tcWaiX'). 

We can similarly evaluate the Gaussian terms in 
the rep resentation of Vps- Defining the integrand in 



Eq. |D1| , ^{u) = {(f>i(f)j\exp{—u^ri2)\4>k4>i)j we see that 
EPL{tc.)Wa = E^a-^K), where u^/p = tl/{l - 
tl) and = i2/^)[{p + u^)^/yp]eMXuyip + 

u'^)]Wa{X). One may directly modify this finite-point 
quadrature summation to evaluate two-electron integrals 
for the Gaussian terms, replacing n with nc, with q, 
and with di. We have applied these modifications to 
the two-electron integrals within the ARGOS code, part 
of the COLUMBUS program package (sT). 



9 



[1] B. Klahri and J.D. Morgan III, J. Chem. Phys. 81, 410 

(1984). 

[2] R. N. Hill, J. Chem. Phys. 83, 1173 (1985). 

[3] J. D. Morgan, in Numerical Determination of the Elec- 
tronic Structure of Atoms, Diatomic, and Polyatomic 
Molecules, edited by M. Defranceschi, and J. Dclhalle 
(Kluwer Academic Publishers, Dordrecht, 1989) p. 49. 

[4] W. Kutzelnigg and W. Klopper, J. Chem. Phys. 94, 1985 
(1991). 

[5] W. Kutzelnigg and J. D. Morgan III, J. Chem. Phys. 96, 

4484, (1992). 

[6] M. Hoffmann-Ostenhoff, T. Hoffman-Ostenhof, and H. 

Stremnitzer, Phys. Rev. Lett. 68, 3857 (1992). 
[7] T. Helgaker, W. Klopper, W. Koch, and J. Noga, J. 

Chem. Phys. 106, 9639 (1997). 
[8] T. Loiningor, H. StoU, H-J. Werner, and A. Savin, Chem. 

Phys. Lett. 275, 151 (1997). 
[9] M. Nooijen and R. J. Bartlett, J. Chem. Phys. 109, 8232 

(1998). 

[10] I. Silancs, J. M. Ugalde, and R. J. Boyd, J. Mol. Struct. 

(THEOCHEM), 527 27 (2000). 
[11] T. Kato, Commun. Pure Applied Math. 10, 151 (1957). 
[12] E. A. Hylleraas, Z. Phys. 51, 347 (1929). 
[13] R. Bukowski, B. Jcziorski, and K. Szalewiczi, J. Chem. 

Phys. 100, 1366, (1994). 
[14] J. D. Baker, D. E. Fround, R. N. Hill, and J. D. Morgan 

III, Phys. Rev. A 41, 1247 (1990). 
[15] T. L. Gilbert, Rev. Mod. Phys. 35, 491 (1963). 
[16] H. F. King, Thcor. Chim. Acta 94, 345 (1996). 
[17] D. R. Hamann, M. Scliliitcr, and C. Chiang, Phys. Rev. 

Lett. 43, 1494 (1979). 
[18] W. E. Pickett, Comp. Phys. Rep. 9, 115 (1989). 
[19] D. R. Hamann, Phys. Rev. B 40, 2980 (1989). 
[20] A. Savin, in Recent Developments and Applications of 

Modern Density Functional Theory edited by J. M. Sem- 

inario (Elsevier, 1996) p. 327. 
[21] C. Filippi and C. J. Umrigar, J. Chem. Phys. 105, 213 

(1996). 

[22] C. J. Umrigar, K. G. Wilson, and J. W. Wilkins, Phys. 

Rev. Lett. 60, 1719 (1987). 
[23] J. C. Greer, J. Chem. Phys. 103, 1821 (1995). 

[24] ,1. C. Grccr, ,1. Comp. Phys. 181, 146 (1998). 
[25] M. Dupuis, J. Rys, and H. F. King, ,]. Chem. Phys. 65, 
111 (1976). 

[26] E. R. Davidson, S. A. Hagstrom, S. J. Chakrvorty, V. M. 

Umar, and C. F. Fischer, Phys. Rev. A 44, 7071 (1991). 
[27] T. H. Dunning Jr., J. Chem. Phys. 90, 1007 (1989) 
[28] A. Schaefcr, H. Horn, and R. Ahlrichs, J. Chem. Phys. 

97, 2571 (1992). 
[29] R. Krishnan, S. J. Binkley, R. Sceger, and J. A. Pople, 

J. Chem. Phys. 72, 650 (1980). 
[30] P. O. Widmark, P. A. Malmqvist, and B. O. Roos, Theor. 

Chim. Acta 77, 291 (1990). 
[31] R. Shcpard, I. Shavitt, R. M. Pitzer, D. C. Comeau, 

M. Pepper, H. Lischka, P. G. Szalay, R. Ahhichs, F. B. 

Brown, and J. Zhao, Int. J. Quant. Chem. Symp. 22, 149 

(1988). 



10 



